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Abstract 

We calculate the emissivity for the pair production of helicity-flipped neutrinos, 
in a way that can be used in supernova calculations. We also present some simple 
estimates which show that such process can act as an efficient energy-loss mechanism 
in the shocked supernova core, and we use this fact to extract neutrino mass limits 
from 5N1987A neutrino observations. 
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1 Introduction 


The observation of neutrinos from SN1987A [1,2], in fair agreement with predictions 
from supernova models, has been used by several authors to bound the properties and 
interactions of various exotic and non- exotic particles [3, 4, 5, 6, 7, 8]. For those particles 
which act as an efficient energy loss mechanism for the supernova, the simplest constraint 
on their masses and couplings can be derived from the fact that the total energy carried 
away by them cannot be greater than the available energy of the star. The detected 
neutrino flux has substantiated several features of supernova theory which are now ac- 
cepted as being standard [9,10]. First, the total emitted energy is 2 — 4 x 10 53 ergs . 
Secondly, it is emitted in the form of neutrino- antineutrino pairs of all species (with 
roughly equal amounts carried by each) formed inside the core via the Z° exchange 1 
process e + + e~ — * v -f P. Finally, these neutrinos are trapped in the core and undergo 
slow thermal diffusion for several seconds until they reach the neutrinosphere, where they 
are released in large quantities. 

If neutrinos are massive Dirac particles, then it is possible to produce neutrino- an- 
tineutrino pairs via the above process such that one of them is non-interacting, i.e either 
a positive helicity neutrino or a negative helicity antineutrino, which does not undergo 
diffussion but leaves the core much faster than its trapped partner 2 . As will be shown, this 
process can lead to significant energy loss on a much shorter timescale from the shocked 
core, depending on how massive the neutrinos are. Our main aim in this paper is to pro- 
vide suitable expressions for the emissivity of this process in a form which can be easily 
incorporated into realistic supernova models to evaluate the energy lost in the form of 
these flipped neutrinos. To illustrate that the process can have significant consequences, 

we also derive, using general considerations, approximate expressions for these energy 

1 Electron type neutrinos are also produced via the W exchange channel. 

J We confine our discussion to the Standard Model, [11,12,13] minimally extended to include Dirac 
masses for the neutrinos. A right-handed Majorana neutrino interacts in a manner similar (but not 
identical) to a right-handed Dirac anti-neutrino and hence cannot provide an avenue for rapid cooling. 
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losses. The emissivity of the process is proportional to the mass squared of the neutrino 
and the seventh power of the core temprature. We show that, using these expressions 
, it is possible to reliably exclude neutrino masses in the range 1 MeV — 100 MeV . The 
above range of 1 — 100 MeV is obtained using low core temperatures [21], and hence is 
on the pessimistic side. If higher temperatures are used, as is the case in some supernova 
models [22] then it is possible to exclude all neutrino masses between 100 KeV - 100 MeV 
using the considerations discussed below. As we will discuss, these mass limits would be 
improved by performing full supernova calculations. 

This paper is organized as follows. In Section 2, we analyze the production amplitudes 
for helicity- flipped neutrinos. In Section 3, we calculate the emissivities associated to 
this processes, in a way that can be used in supernova calculations, and we make some 
estimations about mass limits. 3 . We end by summarizing our most relevant conclusions. 

2 Pair production of helicity-flipped neutrinos 

In this section we obtain the cross-section for the Z° exchange process 

e+ (?i) + e (Pa) — ♦ Ai) + Aj) (2.1) 

Here the ki and pi are the particle momenta and A is the neutrino helicity. We recall that 
within the context of the standard theory, a neutrino interaction eigenstate (vl or vr) 
is a superposition of helicity (A) eigenstates v±, where A = <r.p = ±1. For a relativistic 
particle, this translates into the statement that a vl is predominantly in the A = — 1 state 
and a vr is predominantly in the A = -fl state, with small admixtures of the opposite 
helicity, of order mj E u . In particular, we are interested in the case where the final state 
neutrino has positive helicity, and is ’almost’ non-interacting. 

3 Current experimental bounds on neutrino masses from particle accelerators are < 18eV, m„ < 
0.25MeV, m Vr < 35 MeV [14]. More stringent bounds on the masses of stable neutrinos have been derived 
from big bang cosmology [15,16,17]. The neutrinos we consider here would reasonably be expected to 
decay via mixing and other modes. Bounds from cosmology and astrophysics on the masses and lifetimes 
of unstable neutrinos have been derived in [18,19,20]. 
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In what follows below, the subscipts 1 and 2 denote antiparticle and particle respec- 
tively. The amplitude for ( 2.1) is given by: 

A = ^|fi(fc a , A 2 )7' i (l - 7 5 )v(fc!, A x ) • r(pi)7^(Cv - C A y s )u(p 2 ) (2.2) 

The u and v are the usual Dirac spinors, and use has been made of the fact that, when 
the processes (2.1) occur in the core of the collapsed star, the center of mass energies are 
at most 1 GeV , hence the amplitude may be written in its effective four-fermion form. 
The electron and positron helicity indexes have been supressed since they will be averaged 
over. We then write 


\\\ A \ 2 = (2 - 3) 

(where the spin averaging factors have been explicitly shown) with 

JV *' = J TtKM + m)(l + 7 s - 7 S )(M - rn)(l + 7 ‘ Ab"(l - 7 S )1 (2.4) 

and 

E*" = Tt[(& + M)-r“(Cv - C t 7 s )(b - M)Y(Cv - C J7 S )] (2-5) 

Here &i and s 2 are the spin four-vectors associated with the anti-neutrino and neutrino 
respectively, while m and M are the neutrino and electron masses. These spin vectors 
satisfy the Lorentz invariant conditions 

— 1 j 3, • ki — Oj (2.6) 

and for a relativistic neutrino the additional constraint 

i*i||Aifc t - for t = 1,2 (2.7) 

holds, where 

k* = ( E v , \k(k) (2.8) 
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with k being a unit vector along the three-momentum of the neutrino . 

We now introduce two four-vectors associated with the neutrino pair: 

K? = *5* + msl ; K£ = k$ - ms»- (2.9) 

In conjunction with the properties in ( 2.6) and ( 2.7), these will allow us to write the am- 
plitude squared for the process under consideration in a compact and physically revealing 
form. As a first step towards this, we note that the spin vector may be expressed as 



= 16 (C v + C A ) 2 (p 1 .K 1 )(p 2 -K 2 ) + 

16(CV - C A ) 7 { Pl • K 2 )(p 2 • Ki) + 

16 {Cl - CDM 2 ^ • K 2 ) (2.14) 

Here Cy, Ca are the usual weak vertex factors. From this expression and equations 
( 2.11) and ( 2.13) above one sees that the amplitude vanishes for massless neutrinos , 
as it should . Further, the expression ( 2.14) is akin to the usual weak pair production 
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amplitude with the replacement K; — ► Finally, the flip and non-flip cross-sections are 

related by more than just a simple factor of m 2 /AEl , since the 77 , carry a sign affixed 
to the three-momentum, which in general depends on the nature of the final state (i.e 
whether it is a particle or anti-particle) and its helicity eigenvalue [5] . 

In the next section we proceed to evaluate the emissivity for this process using ( 2.14). 

3 Calculation of emissivity 

We now proceed to calculate the emissivity associated with the pair production of helicity- 
flipped neutrinos. This is given by 

<?m = 4 / ^ (3-1) 

In the last equation Q m is the emissivity, p x (p 2 ) are the e~ (e + ) momenta, E e - and E e + 
are their energies, and / e - (/ e +) are their respective equilibrium Fermi-Dirac distribution 
functions, which are assumed to be 

/e " (£e ' } = l + exp(£ e -- M e-)/r 

fe ^ Ee+ ^ = l + exp(£ e+ + fi,-)/T 
Here, p, e - is the chemical potential for the electrons and T is the temperature (we take 
Kb = 1 for the Boltzmann constant) 

In eq. (3.1), e(pi, P 2 ) stands for 

<Pi,Pi) = J k °2 v MP1P2 -*■ fci fca) (3.4) 

where v is the relative e + — e~ velocity, and da(pip 2 —* k x k 2 ) is the differential cross 
section for the process described in ( 2 . 1 ) (later on, we will also include the process e + e _ — ► 


(3.2) 
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ul^Li the cross section being the same that for e + e~ — > vrvr). Since the typical energies 
which are involved are in the order of magnitude of ~ lOOMeF, electrons will be treated 
as relativistic. Therefore, we have 


Q m = J dx 2 dx 2 x\x\ f t -(xi)f e +(x 3 ) Jocose) 


t(Pl,P2) 


(3.5) 


Where x x = E e - /T, x 2 = E e + /T and 6 is the angle formed by p x and p 2 . 

We have performed the integrals in (3.4). After a lenghty (although strightforward) 
calculation, we obtain 


e(/i>P2) = tyPuPi) (3.6) 

The function ft{px,p 2 ) is given by 

n{p u p 2 ) = E {{Cl + C\) [C x (a) + C a (a)Z^( 1 - cos 6f + C 3 (a)(l - cos 0)] + 

2C v C A C 4 (a) Ee ~ - e E ' + (1 - cos *)} (3.7) 

where E = E e - + E e + , a = |p x + p 2 \ /E and 


. . [(a 6 — a 4 — 4a 3 + 3a 2 4- 4a — 3)//(a) — 2a 5 4- 8a 4 — 4a 3 — 8a 2 4- 6a] . . 

C x (a) = — t 3 - 8 ) 


. [(3a 4 — 4a 3 + 2a 2 — 4a + 3)//(a) — 6a 3 4- 8a 2 — 6a] 

= 2a* 


(3.9) 


^(a) — 


[(a 6 — 3a 4 + a 2 + 4a — 3 )/j(a) + 2a 5 4- 4a 4 — 8a 2 + 6a] 

— 


(3.10) 


C 4 (a) = 


[(3a 4 — 4a 3 4- 2a 2 — 4a 4- 3)/j(a) — 4a 5 — 6a 3 + 8a 2 — 6a] 


(3.11) 
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In the last equations, //(a) = log(^). The integral in cos 6 can be performed numer- 
ically, for different values of E e - and E e +. We realized that, for our purposes, the result 
of this integral can be approximated by the simple expression 


J d { cos 6) 




= Ai (Cy 4- C\){x i + X 2 ) — \iCvCji(x 1 — X 2 ) 


(3.12) 


Here, Ai and A 2 are two numerical constants to be determined later. By inserting (3.6) 
and (3.12) in eq. (3.5) one can get a final expression in terms of the relativistic Fermi 
functions F n (r}), defined as 


. . 1 x n dx 

FM = n\h 1 + exp(® - 77 ) 

here, 77 is the electron degenerancy parameter (77 = n e - /T). We compared the resulting 

expression with the exact numerical integral, as given in (3.5). We then found that 

Ai ~ 2.5 and A 2 — 8.0 is a good choice for 77 < 10, whereas Ai ~ 3.5 and A 2 ~ 9.0 seems 

more appropiate when 77 > 10 . In that way, we obtain the final expressions 





with the following approximations for #( 77 ) : 


( 3 . 14 ) 


H(rj) = 2.5 (Cl + C 2 a )A + ( 77 ) - 8C v C a A-(t,) (3.15) 


(77 < 10 ) 


H( 77 ) = 3.5 (Cf -f C 2 a )A + ( 77 ) - 9C v C A A-( n ) (3.16) 


(77 > 10 ) 
where 
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•4 i (<?) = fifoM-ij) ± F^n)F,{-v) (3.17) 

In eq. (3.14) we introduced a factor of 2 to account for the combined emissivity of the 
two processes e + e“ — ► vlvl and e + e“ — ► vrvr 

Substituing for the values of the constants, one obtains 

Qz- = 5.281 x 10 29 (m/lMeVY(T/l0MeV) 7 H(T]) (3.18) 

erg.cm- 3 s- 1 

In Fig. 1 we show the two approximations given in (3.15) and (3.16) (dashed line and 
dot-dashed line, respectively) for We also have plotted the values that result from 

a direct numerical integration of (3.5). The error in using this approximate expressions is 
only a few percent. 

As one can see from eq. (3.14-17) the emissivity is highly dependent on the tem- 
perature and electron degenerancy in the supernova core. Therefore, in order to get an 
appropiate neutrino mass limit, one should incorporate our expressions in a realistic su- 
pernova calculation. However, in order to motivate such a calculation, we present some 
simple estimates. The highest temperatures in a supernova collapse are reached in the 
shocked, outer core, during the first few seconds of the neutron star cooling. In this re- 
gion, the temperature amounts to several tens of MeV, and r) ss 0 (see ref [21]. For this 
model, the peak temperature is higher than 35 MeV). 

In Fig. 2, we show our estimates for the total luminosity of the shocked core in the 
form of heli city-flipped neutrinos, for various temperatures (we assumed a characteristic 
core radius R = 10 Km). As can be seen, for modest temperatures T = 35 MeV [21] 
one can exclude neutrino masses larger than about 1 MeV , using arguments identical to 
those in ref. 5. However, for supernova models with higher temperatures T ~ 70 MeV 
[22,23], the corresponding mass limit is about lOOA’eV’ (we are convinced that this limits 
can be improved by performing realistic calculations, since the neutrino luminosity would 
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be integrated over a few seconds). 

We now address the question of whether very massive neutrinos can be ruled out 
from the above considerations. For these neutrinos, the re-flipping process into standard 
neutrinos can proceed via scattering with targets such as electrons, neutrons and protons 
in the core [5]. In fact, for neutrino masses larger than a few MeV, the mean free path 
associated with re-flipping becomes comparable to the core radius, and hence one could 
claim that flipped neutrinos will be trapped, rather than freely escaping. However, they 
will be trapped only temporarily, because neutrinos will continue to be flipped. The 
important criterion is when transport is more effectively done by flipped neutrinos rather 
than unflipped ones. Hence, the relevant quantity is the difussion time-scale associated 
with re-flipping 


3 K 2 n<r fu p 

*<*»// — Z2Z 

ir‘c 


(3.19) 


with n being the number density of the targets, and <Tfu p the cross section for re-flipping 
(which is proportional to the neutrino square mass). If tdiff is of the order of 1 sec. or 
so, helicity-flipped neutrinos become effectively trapped, and no longer act as an energy 
loss source. For standard values in (3.19), this only happens when the neutrino mass is 
m > 100 MeV [5]. Therefore, neutrinos which are more massive than this limit can not 
be excluded by helicity-flipping processes (however we note that such high masses for /x 
and r neutrinos are anyways ruled out by accelerator limits and hence reflipping is not a 
relevant concern here). 


4 Conclusions 


We have calculated accurate expressions for the emissivity due to pair production of 
helicity-flipped neutrinos with a Dirac mass. We showed that this process can act as an 
efficient energy-loss mechanism in the shocked core of a supernova. Therefore, this can 
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be used to put limits on the neutrino mass, by using the data of the detected neutrinos 
from SN1987A. 

Because of the high temperature-dependence of the corresponding emissivity, these 
mass limits are better in the case of supernova models with large core temperatures 
[22,23]. If these models are reliable, all Dirac neutrino masses in the range lOOA’eF- 
lOOMeV can be ruled out by our simple estimates. A full supernova calculation should 
certanly lower the 100 KeV limit (probably within a factor of 3 or so), since the neutrino 
luminosity would be integrated over several seconds. In this way, mass limits which are 
comparable to the ones obtained from scattering helicity- flipping processes [5,7] can be 
reproduced by an independent mechanism. 

If one addopts models with much lower core temperatures [21], the corresponding 
(pessimistic) mass limit is higher in about a factor of 10. However, even in this case, the 
extreme sensitivity of the supernova explosion to neutrino properties (specially for the 
so-called delayed explosion scenario [22,24]) could lead to a much more stringent limit. 
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Figure Caption 

Fig. 1 .- The function H(tj) (solid line), as given by a numerical integration of (3.5). 
Also shown, the two approximations corresponding to eq. (3.15) (dashed line) and (3.16) 
(dot-dashed line). 

Fig. 2.- Total emissivity for helicity- flipped neutrinos in the core, as a function of the 
neutrino mass, for three different core temperatures. A typical core radius R = 10 Km 
has been assumed. 
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